INTELSAT IV IN-ORBIT LIQUID SLOSH TESTS AND PROBLEMS 
IN THE THEORETICAL ANALYSIS OF THE DATA 

Victor J. Slabinski 


ABSTRACT 

Each INTELSAT IV dual-spin spacecraft carries liquid hydrazine in four coni- 
spheric propellant tanks mounted on the spinning rotor section. The destabil- 
izing effect of the liquid on attitude-nutation stability was determined from 
an extensive series of in-orbit tests. The liquid slosh driving frequency ratio 
(rotor nutation frequency /rotor spin rate) was varied over the range of 0.58 to 
0.70 for the tests by rotating the spacecraft antenna platform at different rates 
in inertial space. A rotor-mounted accelerometer sensed the spacecraft nutation. 
The observed time constant for the nutation angle increase or decrease was cor- 
rected for the stabilizing contribution of the platform-mounted pendulum damp- 
ers to yield the net destabilizing dedamping contribution from the liquid slosh. 

The in-orbit tests show two unexpected maxima in the dedamping contribution 
at driving frequency ratios that vary with the propellant loading. The rotor 
nutation frequency at the maxima is about one-third of the lowest mode liquid 
slosh frequency given by ground test data for unspun tanks, and thus does not 
correspond to a simple resonance of the liquid. Ground tests with spinning sys- 
tems have produced the same maxima, but the phenomenon is not yet understood. 

A preliminary version of this paper appeared in the Proceedings of the CNES- 
ESA Conference on “Attitude Control of Space Vehicles^ Technological and 
Dynamical Problems Associated with the Presence of liquids,” Toulouse, 

France, October 10-12, 1977 (ESA SP-129), pp. 87-102. The final vemion 
appears in the COMSAT Tech. Rev. 8, 1-40 (Spring 1978). 

This paper is based upon work performed at COMSAT-Plaza under the sponsor- 
ship of the International Telecommunications Satellite Organization (INTELSAT). 
Views expressed in this paper are nOt necessarily those of INTELSAT. 


183 



introduction 

Definitions 

A dual-spin spacecraft such as intelsat iv consists of a spinning rotor 
section and a normally despun antenna platform, with the two sections 
connected by a bearing assembly. The bearing axis is the intended spin 
axis for the spacecraft. This is an unstable axis because, by the common 
definition of dual-spin spacecraft, it is not the axis of maximum moment 
of inertia. The bearing axis may exhibit a small motion about the space- 
craft’s rotational angular momentum vector which is the result of two con- 
ing motions, a steady-state motion at a cone half-angle dw called the wobble 
angle, and a transient motion called nutation whose instantaneous cone 
half-angle 6n is the nutation angle. 

The INTELSAT iv spacecraft carry liquid hydrazine in four partially filled 
propellant tanks mounted on the rotor section. The liquid motion in 
response to nutation tends to increase the nutation angle at the rate kiQx, 
while the platform-mounted pendulum dampers tend to decrease the 
nutation angle at the rate kodN- small nutation angles, the nutation 
angle is observed to vary exponentially with time t\ 

Bn = exp 1^— “7^] (1) 

where t is the time constant for nutation angle decay (t positive) or growth 
(t negative). Its inverse can be expressed as the difference of two positive 
numbers : 

~ = ka- kL ( 2 ) 

T 

where ki, is the damping contribution from the platform dampers and ki 
is the dedamping contribution from liquid slosh on the rotor. The latter 
contribution will be called “dedamping” for brevity. Spacecraft attitude 
stability requires ki < kn so that r is positive and any initial nutation 
damps out. A bound on ki for ail expected operating conditions is thus 
very important for successful spacecraft design and operation. 
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The liquid motion in the tanks is commonly referred to as “fuel slosh.” 
This term is not strictly correct here, since hydrazine is a propellant rather 
than a fuel. Fuel is matter that produces power by combustion, that is, 
by combination with oxygen, while hydrazine produces its thrust on 
INTELSAT IV by Catalytic decomposition. For this reason, the hydrazine 
motion will be called liquid slosh. 

An important determining parameter for liquid slosh effects and hence 
ki is the driving frequency ratio, p, defined as 





where fd is the liquid driving frequency due to nutation, and Wrz is the rotor 
spin rate about the bearing axis, both expressed in units of rad/s. The 
importance of this ratio may be seen by considering that the centrifugal 
acceleration acting on the liquid in the tanks is (4zR, where R is some 
mean liquid distance from the spin axis. The restoring force on the liquid 
when displaced from equilibrium will be proportional to this acceleration 
so that the natural frequencies of the liquid are proportional to the square 
root of this quantity, <arzR^^^. Since the resonant frequencies are propor- 
tional to Urz, the ratio of driving frequency to liquid resonant frequency is 
proportional to p. 

A dual-spin spacecraft offers a distinct advantage for the in-orbit de- 
termination of liquid slosh effects because p can be varied by merely ro- 
tating the platform. The liquid driving frequency during nutation is as- 
sumed to be the rotor nutation frequency, K, observed from rotor-fixed 
axes; hence, the variation of p with the platform spin rate <apz in inertial 
space is given by 



where Hz = bearing axis component of the spacecraft’s rotational 
angular momentum (constant to the first order for 
small nutation angles) 

jpzi Irz = platform and rotor moments of inertia about the bear- 
ing axis 

It = spacecraft moment of inertia about a transverse axis 
through the mass center. 

Since a simple spinner can be considered as the rotor of a dual-spin space- 
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craft which lacks a platform, setting Ip^ = 0 in equation (4) gives 


(simple spinner) 


( 5 ) 


which shows that for a simple spinner p is fixed by its mass properties. 


Test motivation 

Prior to the first intelsat iv launch, Hughes Aircraft Co. performed 
ground tests to predict ki for the expected spacecraft operating conditions 
(Reference 1, pp. 15-16). These tests indicated that any nutation would 
quickly decay, and that ki was sensitive to p variations. However, they 
were performed for only a few different values of fill fraction (the per- 
centage of the propellant tank volume occupied by liquid) and did not 
indicate two troublesome properties of ki which were only discovered in 
1974 as a result of orbital operation of the Intelsat iv series of spacecraft. 
These discoveries raised questions concerning the attitude stability of the 
INTELSAT iv-A series of communications satellites then under development. 

The first was that ki is a strong function of the driving frequency ratio, 
p. This was discovered during the loss of antenna despin control by the 
INTELSAT IV F-7 Spacecraft on March 21, 1974, when the antenna platform 
spun up to 25 rpm before recovery; kz decreased twentyfold during a p 
value change of only 0.03. This discovery, that a small change in p can 
increase or decrease kz by an order of magnitude, raised the possibility 
that INTELSAT iv-A, Operating at different p values, might have kz values 
significantly larger than those found previously in inielsat iv orbital 
operation. 

This discovery also indicated a need for detailed knowledge of the kz 
variation over the range of p values possible during an antenna platform 
spin-up. Such data would indicate whether loss of despin control could 
result in a rapid growth of the nutation angle to '~90°, the “flat spin” 
condition from which spacecraft recovery might be difficult. Previous 
liquid slosh tests on the ground had been restricted to the p value corre- 
sponding to the despun platform condition. 

The second discovery was that appeared to become drastically larger 
with increasing fill fraction. The first five intelsat iv spacecraft launched 
had a 70-percent fill in transfer orbit after their initial spin-up. The sixth, 
INTELSAT IV f-8, was launched on November 21, 1974, with a 76-percent 
fill. The extra propellant was carried as ballast to compensate for an over- 
sized solid propellant apogee motor. Nutation tests in transfer orbit 
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showed that this extra 6-percent fill doubled Ul. Since this fill fraction 
region had not been included in the ground tests, this raised a question 
concerning the magnitude of ki for the intelsat iv-a spacecraft which 
were to use the same propellant tanks with an initial fill of 82 percent. 

To resolve these questions, Comsat conducted an extensive series of 
in-orbit nutation tests on the nearly identical spacecraft in the intelsat iv 
series in spring 1975. These tests took place in a 0-g environment and 
thus avoided the uncertainty inherent in ground testing regarding the 
correction, if any, which is required for the 1-g field. 

Notation test deserlptlon 

The basic test ptocedure consists of inducing spacecraft nutation and 
recording the continuously telemetered output of a nutation accelerometer 
mounted on the rotor. Figure 1 shows a typical nutation trace for a despun 
platform. In this case, nutation has been induced by firing an axial thruster 
for 0.5 s. The sensed acceleration from the nutation oscillates with the 
rotor nutation frequency, Xr. This frequency may be determined directly 
from the time code on the strip chart by noting the time interval for a 
given number of nutation cycles. Scaling the amplitude at selected points 
along the nutation trace and fitting equation (1) to the measurements 
gives the time constant t for the nutation decay or growth. The computed 
platform damper contribution kn is then used in equation (2) to obtain 
the liquid slosh dedamping itz,. 

The in-orbit tests experimentally determined kz as a function of driving 
frequency ratio, p, and propellant tank fill fraction. The driving frequency 
ratio was varied by spinning the antenna platform at different rates in 
inertial space over a range of ±20 rpm. Obviously, the spacecraft did not 
carry communications traffic during the tests. Spacecraft with greatly 
different fill fractions were available, since each satellite of the intelsat iv 
series had been in orbit for a different length of time and had depleted its 
maneuvering propellant by a different amount. Also, intelsat iv f-8 was 
tested before and after large inclination maneuvers to obtain data at fill 
fractions of about 66, 61, and 56 percent from the same spacecraft. Finally, 
to test the initial intelsat iv-a fill fraction, the last spacecraft to be 
launched in the series, intelsat iv f-1, was put into orbit loaded with an 
extra 22 kg of hydrazine propellant to give an initial fill fraction of 82 
percent. Most of this extra propellant was expended in transfer orbit to 
help place this heavily loaded spacecraft into geostationary orbit. 

The test data obtained and attempts to understand it are the subjects 
of this paper. The in-orbit test results prompted an extensive series of 
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Figure 1. Nutation Accelerometer Output for a Despun Platform 




liquid slosh ground tests by the spacecraft manufacturer, Hughes Aircraft 
Co, The large dedamping contributions found in both series of tests 
resulted in important modifications to the intelsat iv-a spacecraft design. 

Basie nutation theory 

The nutation theory presented herein allows for a cross product of 
inertia about the platform bearing axis. This cross product introduces 
additional acceleration frequencies in the rotor propellant tanks. Such 
additional frequencies due to spacecraft asymmetries may introduce un- 
suspected slosh effects. 

Spacecraft model and platform coordinates 

For the present analysis it may be assumed that all damper masses are 
caged and that the propellant is frozen in its equilibrium position so that 
each spacecraft section acts as a rigid body. It may also be assumed that 
the spacecraft is statically balanced so that the mass center of each section 
lies on the bearing axis. The spacecraft mass center then lies at a point O 
on the bearing axis which is taken as the origin of a right-handed Cartesian 
coordinate system XpYpZp fixed in the platform. The Zp axis lies along 
the bearing axis, and the transverse axes XpYp are chosen so that there is 
no platform cross product of inertia about the YpZp axes. The tensors 
giving the moments of inertia of the platform and rotor, respectively, are 



® ^pxz 


(« = 


0 


ipt 


0 


( 6 ) 


Fih 0 




(/r) 


0 Irt 0 


(7) 



0 



where the p subscript outside tensor or vector brackets indicates platform- 
fixed axes. Equation (6) allows for a platform cross product of inertia 
since this cross product produces additional acceleration frequencies 
in the rotor. Equation (7) assumes that the rotor is dynamically balanced 
so that the bearing axis is a principal axis. Both equations assume equal 
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moments of inertia Ij>t and ht about any transverse axis through O so that 
the spacecraft moment of inertia about any transverse axis through O is 
given by 

It = /pt + Irt (8) 

regardless of the relative orientation between platform and rotor. Another 
useful quantity is the spacecraft moment of inertia 4 about the bearing 
axis, 

h — Ipz + Irz (9) 

With respect to the platform-fixed axes, the angular velocities of the 
platform and rotor sections, respectively, in inertial space are given by 


Wp — 


CO 


px 


CO 


PV 


L. ^pz J p 


( 10 ) 


COy 


CO 


pz 


CO 


'py 


0) 


rz I p 


( 11 ) 


since the bearing constrains both sections to have the same Xp and Yp 
components of angular velocity, cop^ and copp. The +Zp direction is chosen 
so that (jirz > 0. The angular momentum, H, of the spacecraft about its 
mass center O is then 


i/ = 


Hz 


H. 


I— _J p 


= (Ip) Wp + (7r) W, 


( 12 ) 


H = 


ItUipz -T Jpxz^pz 

JtO)py 

Jpxz^px ^ Jpz^pz Irz^f 


(13) 
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Because the platform-fixed axes rotate with an angular velocity wp 
relative to inertial axes, the time derivative of any vector such as H is 
given by 


4E = f dH 

dt L dt 

where t is the time, and the brackets around the first term on the right side 
of the equation indicate a time derivative with respect to rotating axes. 
In the absence of external torques, dH/dt — 0 and equation (14) gives 

[#]/=-“' (15) 

as the differential equation of motion. In subsequent equations it will be 
assumed that the despin motor between the platform and rotor sections 
maintains a constant relative spin rate ^ between platform and rotor: 

^ = &Jp8 — oirz (16) 

so that 


+ cop X H 


(14) 


Wp* — <0r» (1*^) 

in equation (15), where a dot over a symbol indicates a time derivative. 

Steady spin in platform coordinates 

The spacecraft model described herein allows steady spin in platform 
coordinates as given by the following particular solution of equations 
(13) and (15) : 

(18) 
V 

o)rz — <^rzd = constant (19) 

where ^pzo(Jf t^pxo 

~ ^px^iJpxz^pxQ Jpz^pzQ “ 1 “ Irz^rz^ 0 , ( 20 ) 


COpO = 


0 

L _1 
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This equation has the solution 


WpiO _ — 1 rt\/l + {b(>}pzo/(>}rzOy 

0)rz0 b 

where h is a dimensionless parameter defined as 


( 21 ) 


“^Ipxif^rz^ ( 22 ) 

tr^rzO “{” (-^pz -^t) ^pzO 

This paper is concerned with spacecraft for which \Ipxz/Irz\ 1 so 
that usually |^| <$C 1, Then, the ambiguous sign in equation (21) may be 
taken as positive to obtain a spin vector nearly parallel to the bearing 
axis. (The negative sign gives the “flat spin” solution.) In this case a 
binomial expansion of the square root in equation (21) gives 


bat 


«p*0 — 




2co„o 


(23) 


COpasO 



Irs^rzO “i" ijpz ^pzQ 


(24) 


Equation (24) shows that, if the platform is spinning (wp^o ^ 0), a nonzero 
platform cross product of inertia gives apx 0 so that steady spin of the 
spacecraft occurs about a platform-fixed axis that is inclined to the bear- 
ing axis. The corresponding rotational angular momentum vector is also 
fixed in the XpZp plane at an angle dw from the bearing axis as obtained 
from equations (12), (13), (18), and (19): 



tan dw ' 

V»1 + H; 
l«..| 


(25) 


tan dw = 

ItQ)pxO Ipxgff^pzQ 


(26) 


Hz, 


where 

JfzQ — Ipxz^pxQ Ipz^pzQ “1" Jrzf^rzO • 

(27) 


Since H is fixed in inertial space, the bearing axis in this motion describes 
a cone of half-angle 6w about ii in inertial space. This is th>- wobble motion. 
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Equation (26) may be expressed in a simple form when }6j <K 1 so that 
equation (24) holds. Equation (24) gives 


fpx^PxO 




r 


0)pzQ L 


IrzODrzQ ”1” ijpz -^0 


] 


(28) 


which may be used to replace the second term in the numerator of equation 
(26) to yield 


tan 6w — 


COpz 0 V Hi 0 



(29) 


Equation (27) shows that the expression within the braces is very nearly 
1 for a small platform cross product ; hence, equations (29) and (24) give 


tan $w ^ 


0)pzQ 


(30) 


tan 6w 


Ipx^pzQ 

Ir^rzO H“ (7p« It) OipzO 


(31) 


as useful forms for numerical computation. Since Hy is zero here, equa- 
tions (25) and (30) show that wp and H are parallel as required by equation 
(15) for steady spin. 

Natation in platform coordinates 

For a dual spin spacecraft, \Ipxz\/It and \Ipxz\/h are quantities of the 
first order of smallness. To sufficient accuracy for this paper, the nutational 
motion about the steady spin motion is given by 


U)p’ = 


^pxo t Ci)t eos Xp 


W( sin Xj 


1^ COpzO ^iO)i COS Xp 


(32) 


with 


Wr« COS Xp 



Ipxz Xp "1“ (^pzQ 

Iz Xp 


(33) 

(34) 
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and 


(35) 


Xp *“ 


Xpo + Xp(? — to) 


where Xp is the platform nutation frequency given by 


• HzO 

Ap y (i)pzQ 


(36) 


Xpo and 03i are constants of integration, and Xp essentially gives the azimuth 
of wp about the steady spin solution for platform axes. This solution shows 
that in the presence of nutation both the , platform and rotor spin rates 
exhibit a small oscillation at the platform nutation frequency when 
Ipxz ^ 0, even for the despun platform ease (wpjo = 0). 

The accuracy of equations (32) and (33) may be verified by direct sub- 
stitution into the differential equations of motion, equations (15) and 
(17). These equations give an exact solution for all wt values when Ipxz = 0; 
in this case oipxo = G from equation (24) and co< becomes simply the magni- 
tude of the transverse component of wp. When Ipxz 0, the solution is 
correct for small ut values through terms of the first order in Ito}t/Hzo 
provided that Xp is not small, more specifically, provided that the in- 
equalities 

(37) 


hold for j = 0, 1, 2, 3. The different j values correspond to different 
neglected terms in the differential equations. Physically, the platform spin 
rate ojp^o ^ Hzo/It, which gives Xp Ci:; 0, is excluded because the nutational 
motion assumes a different character which is not of interest for this paper. 

For small wt values, equations (32), (33), and (13) show that H observed 
from platform axes precesses about its steady spin direction at a nearly 
constant angle, the nutation angle % : 


tan ^ 


ItUt 

Hzo 


(38) 


For small angles, is linearly proportional to «t. 

Notatloii in rotor coordinates 

The nutation seen on the rotor is described by using right-handed rotor- 
fixed Cartesian axes XrYrZr with the origin at O and the -\-Zr axis coin- 
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cident with the + Zp and bearing axes. For these new axes, 



COr = 




ry 


L J 


(39) 


where an r subscript outside vector brackets indicates rotor-fixed axes, 
and 

Wfx ^ cos ^ ^pu ^ (^^) 


Ury = Wpx sin ^ + Oipy COS ^ 


(41) 


where ^ gives the angle about the bearing axis between the -\-Xr and -^Xp 
directions. Equations (32), (33), (35), (36), and (17) then give 


where 


(t)r — 


oipxQ COS yp -f- coi COS X 
Wpxo sin yp (at sin X 
co/gfO cos Xp 


Xr Xp + ^ XrO + *Xr(/ 



In equation (43), Xr is the rotor nutation frequency given by 


(42) 

(43) 


Xr = Xp -f ^ 


(44) 


'Xr 


Xr ~ WraO 


where 


Xr = 




(45) 


(46) 


X/ is the inertial nutation rate at sniall nutation angles, the rate at which 
the bearing axis precesses about H in inertial space, ^^ote that terms with 
three different frequencies are present in the rotor angular velocity . 
Equation (45) shows that Xr may be arbitrarily varied in the in-orbit 
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tests. Since Ipxz ~ 0, changing the platform spin rate varies w„o but not 
Xr because Hz stays constant. From equation (27), in the absence of 
thruster pulses, 

H zQ “ ^pz^pzo "4* ^Tz^Tzo ~ J rz^Tzd constant (^*^) 

where Wrzd is the rotor spin rate with the platform despun. Equation (45) 
may be written as 




■- ^rzd 




(48) 


Eliminating aj„o from equations (45) and (47) yields equation (4) for 
the variation of p as a function of Wpjo. An alternate form giving p as a 
function of ^ is obtained by eliminating top*o and o}„o from equations 
(16), (47), and (48): 




(49) 


Rotor accelerations and accelerometer 

The acceleration a about the spacecraft mass center O of a point P 
fixed in the rotor is obtained as follows. The Tr axis may be chosen so that 
the position vector for P is given by 


Rp 


r = 


0 


(50) 


This point has a velocity wr X r about O; hence, for rotating rotor axes, 


(Wr X 7) j + Wr X (oJr X 7) 

(51)* 

r <^T~\ -* /-* -* 2 

1 — J X r + Wr(wr • r ) —r ui . 

(52) 
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Keeping terms to first order in Wpxo and uu in using equation (42) to 
eliminate gives 


[ — oirzO “1“ ^kzfatOirzQ COS Xp] 


' 1 Zp COS ^ I ^r^o) COS 

vpf^x^^Xp sm XpJ 

\ “F Zp lp}pxC0pzO sin Ip -f- Wi(Xr "}■ Oirzo) Sin Xr] 
1 Rp ^0)pxo(/^^rz0 *“ ^pz(^ COS “p 
“1“ — Xr “i“ ^»>rxo) COS Xr] 


(53) 


The large constant term — in the Xr axis component gives the cen- 
tripetal acceleration due to rotor spin. 

Equation (53) makes it possible to interpret the output of an accel- 
erometer located at P. There are two accelerometers mounted on the rotor 
at Rp = 0.99 m with their sensitive axis parallel to the bearing (Zr) axis. 
Filtering their output to pass only frequencies between 0.03 and 35 Hz 
eliminates any contribution from the large Xr axis constant term that 
might be sensed due to accelerometer mis align ment. The filtered output 
is transmitted by fm real-time telemetry to the ground where it is detected 
and displayed versus time on a strip chart recorder. The filter readily 
passes both Zr axis terms for intelsat iv at its typical 50-rpm rotor spin 
rate: the wobble term with frequency ^ ~ — 5 rad/s and amplitude 


z4jp' — Rp(t}pxo(2o}rz0 (*}pzo) 


(54) 


Aw^ 


RpJ PXZ^^^ Tzd Oipzo) 0)pzQ 

Irz^rzO "l~ ijpz ~ WjjzO 


(55) 


from equations (53) and (24), and the nutation term with frequency 
Xr ~ —4 rad/s and amplitude proportional to coj, so that it is a measure 
of the nutation angle. For small platform spin rates, equation (55) gives 
the approximation 

Aw ^ 2Rp ~ uU (56) 

Irz 


which shows, that the wobble term amplitude depends on the ratio Ipxz/Irz- 
Both the wobble and the nutation terms can appear in the nutation 
trace and thereby produce a beat phenomenon when the platform spins 
as shown in Figure 2 for intelsat iv-a, for which Ipzz/Irz = 0.06. This 
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Figure 2. Nutation Accelerometer Output for a Spinning Platform with a 

Large Inertial Cross Product 




moment of inertia ratio is large because of the Intelsat iv-a antenna 
mast design. The resulting wobble term completely dominates the nuta- 
tion trace at high platform spin rates and complicates the hand reduction 
of the nutation strip charts to yield nutation time constants. It is also 
possible that the liquid excitation by a large wobble term results in a 
significant contribution to the observed dedamping, intelsat iv does not 
have these problems because its antenna platform was dynamically 
balanced about the bearing axis to give a small /p*j value. Its observed 
wobble motion gives \Ipxz\ ~ 0.13 kg*rrfi so that Ipxz/Irz ^ 0.0006 for 
this spacecraft. Thus, intelsat iv is ideally suited for spinning platform 
slosh tests. Its wobble term is detectable for [wp^oj > 8 rpm, but becomes 
troublesome only for copzo 18 rpm near the Xp ~ 0 platform resonance. 

Note that some terms in a with frequency Xp ~ 2 rad/s are not detected 
by the rotor accelerometer because they do not appear for the Zr axis. 
This observation indicates that all the rotor acceleration frequencies 
which might produce liquid slosh cannot be determined by merely per- 
forming a spectrum analysis of the telemetry record. 

If P is a point in the propellant tank liquid, the Yr and Zr axis com- 
ponents of a give the driving accelerations parallel to the free liquid 
surface which produce liquid slosh. The terms proportional to wj are 
responsible for the dedamping, since their effect on the liquid increases 
with the nutation angle. Such terms here have two different frequencies, 
Xr and Xp, which are not harmonically related. The Xp terms are not due 
to nonlinear effects, but arise from a theory linear in co«. The Xr terms have 
the largest amplitude by a factor of 100; hence, this paper adopts the 
usual assumption that these terms give the only significant accelerations 
producing liquid slosh. The Xp terms should not be completely dismissed, 
however. Although the Xp term in the Yr component has the small factor 
kz, the liquid response at this frequency could be very large if |Xp] equalled 
a natural frequency of the liquid. Note that this term persists even for a 
despun platform. 

Theoreticians should develop the complete theory of such small accel- 
eration terms due to vehicle inertia asymmetries and mass center mis- 
alignments because of the possibility of frequency coincidence with liquid 
resonances. It is possible that such a frequency coincidence with a well- 
understood liquid resonance might explain the anomalous resonances 
observed in orbit. Similar small acceleration terms should also occur for 
spinning test fixtures on the ground due to interactions with the earth’s 
gravitational field. The amplitudes of such terms do not necessarily become 
zero because of a rapid vehicle spin rate. 
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Nutation test details 
Nutation test procedure 

The INTELSAT IV pendulum dampers are mounted in an open frame 
structure on the antenna platform. As a result, their temperature under- 
goes a large daily variation over a range of 100°C due to changes in 
shadowing and the direction of the incident sunlight as the earth-pointing 
antenna platform rotates at 1 rev/day in inertial space. This temperature 
variation causes a 40-percent change in the eddy-current damping force. 
To eliminate this variation from an extensive series of nutation tests 
extending over many hours, the platform was spun in sunlight in a “ro- 
tisserie mode” at 1 rpm for 6 hours before the tests to establish an equi- 
librium temperature. 

For each test, the antenna despin control electronics were operated in 
the rate memory hold mode (Reference 2, pp. 295-299) and commanded 
to the proper relative spin rate yj/ computed from equation (49) to give 
the desired p value. Spacecraft telemetry sent through an omnidirectional 
antenna indicated when the desired platform spin rate had been estab- 
lished. This spin rate can also be measured directly by recording the signal 
level of a carrier sent through a global beam antenna on the platform; 
the signal strength rises from zero once each platform spin period as the 
global beam illuminates the earth. With the spin rate established, accel- 
erometer data were usually obtained over at least one nutation time con- 
stant. Only for tests with r > 600 s was this not always done. 

For some unstable platform spin rates, the small nutation already present 
from starting transients would quickly increase to measurable levels. For 
other spin rates, measurable nutation was already present from a previous 
test. An axial thruster was fired to induce nutation only when necessary 
to avoid violating thermal constraints on thruster restarts at short time 
intervals. The tests were restricted to nutation angles less than 1°. To 
reduce the nutation angle following tests that gave unstable nutation, the 
platform was returned to the rotisserie mode or to some other platform 
spin rate that gave stable nutation rather than to the despun condition, 
in order to maintain the damper equilibrium temperature. 

For tbe most accurate presentation of the data, the value of Xr used to 
compute the p value actually occurring in each test should be measured 
on the nutation strip chart recording. For most of the tests, Xr was com- 
puted from equation (48) simply because the rotor spin rate is readily 
obtained by telemetry from the rotor-mounted sun and earth sensors. 
Equation (48) gives a useful check on the spacecraft mass properties corn- 
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puted from prelaunch measurements, but it strictly holds only for a rigid 
spacecraft. Liquid slosh in response to nutation may result in measurable 
differejnces from the values of equation (48). These differences could 
provide clues to the nature of the anomalous liquid resonances observed 
in the tests. 

Reference 2 describes the intelsat iv spacecraft in great detail. Table I 


Table 1. Mass Properties for Intelsat IV F-8 



Transfer 

Orbit 

Synchronous Orbit 


Start of 
Life 

End of 
Life 

Spacecraft Mass (kg) 

1,404 

753 

608 

Hydrazine Mass (kg) 

144 

142 

0 

Fill Fraction, f 

0.758 

0.747 

0 

Tank Center Location zr Above 
Spacecraft Mass Center 0 (m) 

+0.24 

-0.24 

-0.39 

Moments of Inertia (kg m*) 





93 

93 

93 

Ir, 

321 

255 

201 

h 

1,131 

702 

663 


lists the mass properties of a typical intelsat iv spacecraft based on ground 
measurements; adequate moment of inertia values for other synchronous 
orbit hydrazine loadings may be found by linear interpolation in the table. 
The following two subsections provide those details on the platform nuta- 
tion dampers and propellant tanks necessary for a detailed understanding 
of the in-orbit tests. 

Platform nutotloii dampers 

Neer [1] adequately describes the two platform-mounted pendulum 
nutation dampers on intelsat iv and gives typical parameter values. The 
damping contribution from the ith damper is given by 


^D.i = 


rtiiZDi 






(57) 
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where m, — mass of the zth damper pendulum 

zd% = distance along the spacecraft bearing axis from the 
spacecraft mass center to the /th pendulum mass 

o)ni = natural frequency of the rth pendulum 

= damping factor expressed as fraction of critical damping. 

Haines and Leondes [3] derive this result using an energy sink model, 
which is an adequate approximation for this paper. The form of the equa- 
tion presented here avoids singularities at Xp = 0. Superposition is as- 
sumed to hold so that 


kj) = Atd.i + ko,2 . (58) 

Equations (57) and (58) were evaluated for each nutation test based on 
preflight ground measurements of the damper parameters. Then the 
liquid dedamping could be evaluated from the observed nutation time 
constant t using equation (2) rewritten in the form 

kL = ko — - . (59) 

r 

In some tests, the dedamping contribution is negligible compared to the 
damping contribution so that the right side of equation (59) involves the 
difference of two nearly equal numbers. The computed ki values are very 
uncertain in such cases; therefore, such values have not been included in 
this paper. 

The pendulum dampers swing in a plane parallel to the bearing axis but 
are not mounted directly on the bearing axis. At platform spin rates with 
absolute values above ~10 rpm, centrifugal forces peg the dampers against 
their stops, which prevents them, from swinging. The most accurate de- 
damping values are determined with the dampers pegged because their 
damping contribution is then zero and ki = — I/t with no error due to 
uncertain damper parameters. The pegged dampers contribute an almost 
negligible 0.2 kg • to Ipxz- 

Conlspherlc propellant tanks and fill fraction 

Four conispheric tanks are mounted on the rotor at intervals of 90° 
about the bearing axis. Two such tanks are shown in cross section in 
Figure 3. The tank interior consists of a spherical section with an internal 
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Figure 3. INTELSAT IV Spacecraft \vith Cutaway View of Rotor Section 
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radius of 0.221 m and a tangential conical section v.'ith a full cone angle 
of 85®. The cone axis passes through the bearing axis at an angle of 45®. 
This conispheric design allows propellant to drain completely from the 
tank through the apex of the cone under two conditions: on the ground 
when the spacecraft is vertical, and in orbit during thruster firings under 
the action of centrifugal forces and nitrogen gas pressurant in the tank. 
In theoretical discussions for fill fractions greater than 40 percent, the 
tanks are often assumed to be spherical since the conical section is then 
nearly full and the free liquid surface lies almost entirely within the spheri- 
cal section. The cone volume outside a sphere inscribed in the tank is only 
3.7 percent of the tank volume. 

There are two independent propulsion sjibsystems; each pair of dia- 
metrically opposite tanks is connected to the same subsystem by two lines. 
A gas pressure equalization line connects at the point on the spherical 
wall directly opposite the cone apex, and a common propellant line con- 
nects at the cone apex so that the two tanks may drain evenly during use 
and maintain rotor balance. There is no oscillating transfer of liquid be- 
tween tanks in response to nutation because of the high flow resistance in 
the connecting propellant line. This resistance is due to the length (3.25 m) 
and small inside diameter (5.3 mm) of the line. A high resistance is con- 
firmed by the very long (18-min) time constant observed for the in-orbit 
liquid redistribution between tanks after an initial imbalance [4]. It will be 
assumed that the two tanks in a subsystem have identical fill fractions. 

The center of each spherical section lies at a distance Rt = 0.578 m from 
the bearing axis, and at a bearing axis position Zf — -1-0.24 m above the 
spacecraft mass center O while in transfer orbit. The expulsion of 650 kg 
of solid propellant during apogee motor firing (AMF) shifts O toward the 
antenna platform so that zt — —0.24 in after AMF. Since the tanks are 
nearly the same distance from O before and after AMF, it should be 
possible to combine dedamping values determined before and after. This 
is desirable because of the different p values obtainable before and after 
AMF; p due to nutation with the platform despun changes from 0.72 
to 0.64. 

The following considerations enter into the computation of the tank 
fill fraction. Hughes Aircraft Co. reports that each intelsat iv tank has 
an internal volume of 0.04689 m® on the average with a standard deviation 
of 0.14 percent. Multiplying this volume by 1.008 X 10® kg/m®, the 
density of hydrazine at the tank operating temperature of 20®C, gives 
47.26 kg for the hydrazine capacity of an unpressurized tank. The initial 
nitrogen gas pressure of 1.7 MN/m® increases the tank volume by 0.10 
percent because of the thin wall construction. This volume increase is 
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reduced as the gas pressure decreases to 0.7 MN/m^ with the depletion 
of propellant. The volume change should be included in the computation 
of the fill fraction for comparison of in-orbit tests with ground test data 
because the ground tests are often made with unpressurized tanks. The 
fill fraction, /, is computed by correcting the propellant mass, nip (in kg), 
in the propulsion system for the propellant in the connecting lines and 
the spin-up thruster propellant reservoir (about 0.3 1 kg for each subsystem) 
to find the mass actually in the tanks; this figure is then divided by the 
total mass that the two tanks can hold at the existing pressure: 


/•= mp - 0.31 kg 

^ ~ 2(47.26 kg) [1 -b (0.0006 mVMN) p\ 


(60) 


where the tank pressure, p (in MN/m^), is given by spacecraft telemetry. 
This formula has been used to compute all the in-orbit fill fractions 
quoted in this paper. 

The propellant mass on board a spacecraft is somewhat uncertain be- 
cause it is based on bookkeeping records of maneuvers and a predicted 
value of propellant flow when a thruster valve is opened. The accuracy 
of this propellant accounting was recently checked at Comsat for some 
INTELSAT m satellites built by TRW, Inc. In May of 1977, the hydrazine 
was used to exhaustion in propelling these satellites into orbits high above 
synchronous altitude where they were abandoned. The end of spin speed 
changes and an accompanying more rapid drop in tank pressure indicated 
the point of exhaustion. Table 2 lists the hydrazine remaining in each 
system at this point according to Comsat’s accounting. 


Table 2. INTELSAT III Hydrazine Use Accounting Errors 


Parameter 

INTELSAT III F-2 
Propellant 
System 
I 2 

INTELSAT HI F-6 
Propellant 
System 
1 2 

Hydrazine at Launch (kg) 

10.91 

10.91 

10.30 

10.30 

Hydrazine Left According to 





Accounting at Exhaustion (kg) 

0.50 

1.13 

0.18 

0.54 

Accounting Underestimate of 





Propellant Use (%) 

4.6 

10.4 

1.8 

5.3 

mean (%) 


5.5 


standard deviation (%) 


±3 

.1 
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Table 2 shows that the accounting underestimated the propellant use 
by an average of 5.5 percent. If a similar error occurs for intelsat iv, 
the fill fractions listed in this paper for the most depleted tanks will be in 
error by 1 percent. 

Nutation test results and interpretation 

Data scaling 

It is first necessary to consider whether any scaling should be applied 
to the observed dedamping values to make the data independent of 
arbitrary test parameters such as the rotor spin rate. For this discussion, 
a particular model for liquid slosh is assumed, not because this model 
accurately reflects reality, but only because it provides a simple illustration 
of scaling problems. Specifically, it is assumed that liquid slosh produces 
its destabilizing effect by the back and forth pendulum motion of the 
liquid mass in response to the nutation, lorillo [5] models the liquid by an 
analogous mass-spring-dashpot damper mounted on the rotor at the tank 
center. This damper will be called the lorillo slosh model. lorillo shows that 
the dedamping from longitudinal oscillation (oscillation parallel to the 
bearing axis) is then given by 


ki 


(Xy fa>r^)(Xr ~ XfCd 

lit (xj - 


(61) 


where Cd is the damping constant (the resistive force per unit velocity), 
rrid is the mass, and Wd is the natural frequency of the damper. A similar 
expression holds for transverse oscillation, with R\. replaced by z\. Be- 
cause Xr is generally negative for a dual spin spacecraft, the driving fre- 
quency ratio is given by 


and 


P 





(62) 


(1 - p)(l -b p)^ pCd 
lit (p^ — Pr-iJ® + [pCd/irndiOrz)]'^ 


where pres = (64) 

gives the resonant frequency ratio. For the liquid case, since wd is propor- 
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tional to «rs as pi^viously explained, pres is independent of Wr-. 

The following special cases give the two extremes for spin rate scaling. 
For small damping constants, theexpression within braces in the denom- 
inator is much smaller than Then, for frequency ratios far below and 
far above resonance. 



j?| (l ~ p)(l H- p)^ pc, 

24 >L 


(P Pres) 


m 



R^ d - p)(l -H p)^ pea 

2/; 


(P Pfe) 


( 66 ) 


and the dedamping is independent of rotor spin rate. This is one scaling 
extreme. At the damj^r resonance. 


= (^ ~ 4- Presf 

2^t Pres^d 


and the dedamping varies as the square of the rotor spin rate. This is the 
opposite scaling extreme. 

These three cases indicate a scaling caveat : the dedamping value meas- 
ured at a single rotor spin rate cannot be sealed to a different spin rate 
without some detailed knowledge of the dedamper mechanism such as 
the relation between the frequency ratio used in the test and resonance. 
Therefore, scaling will not be applied to the observed dedamping values 
with the exception of the following application. Equation (63) indicates 
that the dedamping value kz measured for a spacecraft varies inversely 
with its transverse moment of inertia, It, so that the dedamping k[ for a 
spacecraft with a different transverse inertia // but identical tank geometry 
relative to the mass center is given by 

K = 'i,K ■ m 

This scaling is us^ to combine transfer orbit data with data taken after 

AMF. 

An important scaling application is the spin rate scaling of ground 
measurements using spinning test fixtures to orbit' ng satellites. These 
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ground tests usually involve very high spin rates to obtain a high Froude 
number, Fk (the ratio of centrifugal acceleration to gravitational accelera- 
tion), at the tank in order to minimize gravitational effects, while the 
in-orbit satellite has a much lower spin rate. However, the following diffi- 
culty arises when scaling isolated dedamping values from ground tests 
using full scale tanks. For isolated test points, the relation between the 
driving frequency ratio used and the liquid resonances is not known; 
hence, the scaling caveat just given warns that a simple scaling is not 
possible. In addition, an upper bound on the in-orbit dedamping may be 
obtained by using the scaling extreme from the lorillo slosh model that 
gives the highest (most pessimistic) dedamping value. For scaling to lower 
spin rates, an upper bound is thus obtained by assuming that dedamping 
is independent of spin rate. 

Cautions on the time constant group 

The time constant group is discussed here since this quantity is some- 
times misapplied in spin rate scaling. From a dimensional analysis, Neer 
and Salvatore, in Chapter 4 of Reference 6, propose the dimensionless time 
constant group Teg as a convenient parameter for correlating ground and 
in-orbit dedamping measurements. They give 

TtudlMrztDD 

Teg = function {Re, Fr,/, S, p, dr/Rry Zt/Rt) (70) 

where n = number of tanks 

p = mass density of the liquid 

dr = a characteristic tank dimension such as the sphere diameter 
fdr conispheric tanks 

'’’dd — dedamping “time constant” 

Fr = Froude number at the tank 

= centrifugal acceleration of tanks/gravitational acceleration 
(an infinite ratio in a 0-g field) 

S = tank shape factor 

Rt — distance of the tank (sphere) center from the bearing axis 
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zt — tank height above spacecraft mass center 


ri>x> ^ 7- (71) 

with the Reynolds number Re defined as 

J?. = (72) 

where tj is the liquid viscosity. 

The functional relation in equation (70) is to be discovered empirically 
by experiment. This functional relation is useful only if the dimensionless 
parameters in the argument list can be varied individually. However, for 
the in-orbit dedamping measurements described herein, at least two of 
the parameters are always varied simultaneously ; spinning the antenna 
platform at different inertial spin rates varies not only p but also Re 
through the accompanying rotor spin rate change, as shown by equation 
(72). Therefpre, the present paper does not use the time constant group. 

For applications, equations (69) and (71) yield 


ki 


rtfld^rz 

ItT eg 


(73) 


This equation must always be applied cautiously because Teg can be a 
function of the other parameters on the right side of the equation. For 
example, equation (73) gives no definite theoretical relation between ki 
and u>rz since equations (70) and (72) indicate that Teg is a possible function 
of Urz through the Reynolds number. From ground test data, Neer and 
Salvatore conclude in Chapter 7 of Reference 6 that Teg is independent of 
(j)rz at a resonance (fcz, maximum); then ki varies linearly with the rotor 
spin rate. Equation (67) from the lorillo slosh model predicts a different 
result in this case, that ki varies as the square of the rotor spin rate. If 
the Neer and Salvatore linear dependence could be verified by in-orbit 
tests, it could provide a valuable clue to the nature of the observed de- 
damping resonances. 

Obiserved liquid resonances 

Figures 4-8 show the observed dedamping values plotted versus the 
driving frequency ratio for various fill fractions. Figure 8, which corn- 
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pares three different fill fractions, clearly shows that a change of a few 
percent in the fill fraction can result in a drastic change in the dedamping 
for a given p value. Note especially the increase in dedamping with in- 
creasing fill fraction at p = 0.72, the ratio for a despun platform while in 
transfer orbit. Each figure lists the dates for the nutation tests, the satellite 
used, and the satellite h value at the time of the tests for data scaling appli- 
cations using equation (68). Figure 8 includes transfer orbit data scaled 
to the value // = 700 kg»m- so that they may be readily compared with 
the rest of the data taken after amf. The fill fractions for both propellant 
systems are always given since they may be very unequal. 

PLATR3RM INERTIAL SPIN RATE (RPM) 


+20 +10 0 -10 -20 



Figure 4. Tank Liquid Dedamping Contribution for a 51- Percent 

Fill Fraction 


210 




The data show distinct dedamping maxima which are interpreted as 
liquid resonances. These are called ikis anomalous resonances mvx their 
origin is not known. Two maxima are evident in some of the data: Figure 
4 definitely shows two peaks, while Figure 5 shows one peak and part of 
a second. Figure 9 shows that the frequency ratio pres at a maximum 
moves to higher values with increasing fill fraction, indicating that the 
anomalous resonances are actually due to the presence of liquid hydrazine 
on the rotor. Parallel straight lines have been drawn on the graph to 
indicate the trend of the data. 

Thruster firings during some test series changed the spacecraft rota- 



0.58 0.60 0.62 0.64 0.6 ';i 0.68 0 70 


DRIVING FREQUENCY RATIO P 

Figure 5, Tank Liquid Dedamping Contributi m for a 56- Percent 

Fill Fraction 




tional angular momentum so that different data sets for a given fill frac- 
tion may correspond to a significantly different rotor spin rate with 
the platform despun. Such data sets are distinguished by different plotting 
symbols. Figure 5 compares two such data sets taken before and after a 
rotor spin-up which changed (a„d from 48.54 to 51.74 rpm. The two data 
sets fall along the same curve although the data with the higher o)„d 
value appear to have slightly higher dedamping values. At the resonance 
peak, equation (67) from the lorillo slosh model predicts a 12-percent 
difference between the data sets, while equation (73) from the time con- 

PLATFORM INERTIAL SPIN RATE a>pz (RPM) 
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Figure 6. Tank Liquid Dedamping Contribution for a 61 -Percent 

Fill Fraction 
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stant group predicts only a 6-percent difference. There is too much scatter 
in the data to decide between the two predictions. 

Problems in the interpretation of the 
observed resonances 

Before the in-orbit tests, no liquid resonance was expected within the 
range of frequency ratios tested. This opinion was based on ground 
measurements of the lowest mode natural frequency of the liquid in non- 
spinning spherical and conispheric tanks. Tests at Comsat Laboratories 

PLATFORM INERTIAL SPIN RATE (RPM) 
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Figure 7. Tank Liquid Dedamping Contribution for a 63 67 -Percent 

Fill Fraction 
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Figure 8. Tank Liquid Dedamping Contribution for Three Different Fill 
Fractions Obtained During Launch Operations 


showed that the two tank geometries have the same natural frequency 
within 10 percent for a given fill fraction between 8 and 81 percent. These 
measurements indicated that the liquid resonance frequencies would 
always be much higher than the nutation frequency. 

For small amplitudes, the lowest mode oscillation consists of liquid 
motion in the spherical tank section much like a rigid body undergoing 
a pendulum oseillation about an axis through the sphere center. This 
motion with a nearly planar free liquid surface is similar to the residual 
sloshing motion commonly observed in liquid containers such as water 
buckets and wine glasses. Because this is a pendulum-likp motion, the 


214 





Figure 9. Driving Frequency Ratio at an Anomalous Resonance vs Fill 
Fraction for the Conispheric Tanks 


ground measurernents of natural frequeney collected in Reference 7 can 
be scaled to in-orbit conditions by replacing the acceleration of gravity, g, 
in that figure with the appropriate centrifugal acceleration in the space- 
craft tank. For longitudinal oscillations this acceleration is where 

Rc is the distance from the bearing axis to the liquid mass center, while 
for transverse oscillations the effective acceleration is The latter 

acceleration is smaller because, for transverse displacements from equi- 
librium, the mass center of the equivalent physical pendulum does not 
approach the bearing axis as rapidly as for longitudinal displacements. 
Indeed, for a pendulum pivot on the bearing axis {Rt — 0), transverse 
displacements keep the mass center at a constant distance from the bear- 
ing axis so there would be no centrifugal restoring force to produce oscil- 
lation. Table 3 compares the resulting predictions with the observed res- 
onances for different fill fractions and shows that the predicted values are 
--^3 times the observed values. 

The fact that the predicted frequency ratio for a transverse resonance is 
roughly three times the observed frequency ratio has led to speculation 
that a liquid resonance at a multiple of the nutation frequency was excited 
through nonlinear effects. This superharmonic resonance hypothesis to 
explain the dedamping peaks can be dismissed by the folhnving two argu- 
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Table 3. Observed In-Orbit Liquid Resonances and Predicted 
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♦The radius of the spherical tank section, r = cfr/2 = 0.221 m. The distance from the bearing axis to the tank sphere center, 
V = 0.578 m. 



ments. First, an unpublished analysis by the author shows that, for an 
oscillating mass motion on the spacecraft in response to nutation, only 
the frequency component of the oscillation at the nutation frequency 
contributes to the dedamping. The tank liquid may undergo a large ampli- 
tude oscillation at its natural frequency due to a superharmonic resonance, 
but this oscillation is not at the nutation frequency and thus cannot pro- 
duce a secular change in the nutation angle. Second, a careful analysis of 
several dozen nutation strip chart recordings shows that the observed 
nutation amplitude varies exponentially with time as given by equation 
(1).* Therefore, 

. d 

d = - (74) 

T 

where t is a constant. The time rate of change of the nutation angle is thus 
proportional to the nutation angle; i.e., the damping and dedamping 
mechanisms are operating in their linear range and nonlinear effects are 
not important. 

To illustrate the difficulty of using a liquid pendulum-motion resonance 
to explain the in-orbit data, it can be attempted to model the dedamping 
by a simple pendulum of length Ic pivoted at a point fixed in the rotor with 
the pendulum point mass positioned within the liquid volume. Such 
pendulum models have the desired property that, as in the case of liquids, 
their natural frequency, un, is proportional to the rotor spin rate. For 
longitudinal oscillations, 



where Rc is now the distance of the pendulum mass from the bearing axis. 
For transverse oscillations. 


*Some small departures from an exponential fit to intelsat iv nutation data 
have been observed as follows. The amplitude for the first five cycles following a 
thruster pulse is often less than that predicted by an exponential fit, but this is 
attributed to an initial transient in the damper and liquid motions. The first few 
cycles are therefore excluded from the fit. The residuals i'rom an exponential fit 
to the peaks of the following cycles sometimes show a systematic slow oscilla- 
tion with period ^^30 s, but the amplitude is of the same order as the measure- 
ment accuracy so that this is not conclusive evidence for departure from an 
exponential variation. 
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(76) 



For the ratio Wn/corj to be equal to the observed dedamping resonance 
ratio of '^0.66, longitudinal oseillations require /<; ci- 2.3 /?«; that is, the 
pendulum mass and pendulum pivot must be on opposite sides of the 
bearing axis. This is a physically unrealistic model. Transverse oscilla- 
tions require 4 ~ 0.7 Rc\ that is, the pendulum pivot is outside the tank 
near the bearing axis instead of near the tank center. This is also an un- 
satisfactory model. 

Conclusions and suggestions tor 
turther investigations 

The in-orbit dedamping measurements and the considerations of the 
previous subsection lead to the following conclusions: (a) over the range 
of p values tested, there are two anomalous resonance peaks; (b) their 
cause is unknown. It is highly desirable to discover the excitation mech- 
anism that causes the anomalous resonances so that mathematical models 
can be constructed that would allow ground test results to be easily scaled 
to in-orbit conditions. Useful information could be obtained from a 
spinning test fixture such as the angular motion fixture described by 
Martin in Reference 8 (pp. 9--12) or the spinning air bearing fixture at 
Hughes Aircraft Co. (Reference 6, Chapter 3; also Reference 1, p. 75, 
Figure 7) since these spinning fixtures also exhibit anomalous liquid 
resonances. However, these fixtures must be instrumented to provide 
quantitative measurements of the actual motion of the liquid, A closed 
circuit television camera is sometimes mounted on the rotating test vehicle 
to view the liquid contained in a transparent tank. Such cameras have not 
proven useful in quantitative work because nonlinear effects are likely 
to be introduced by the large nutatidnal motion usually required to give 
measurable liquid motion. Such cameras can give useful qualitative 
results, however, 

The author has proposed th&t liquid height gauges be used to obtain 
quantitative measurements of the amplitude and phase of the rise and 
fall of points on the free liquid surface. Even coarse (10-percent) measure- 
ments could define the surface motion from which the mode of liquid 
oscillation (longitudinal, transverse, or other) at resonance could be 
deduced. The gauge proposed here is based on the electrical conductance 
between an electrical probe and the tap water used in place of hydrazine 
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as the tank liquid in ground tests. As shown in Figure 10, two bare elec- 
trical probes ~2 mm in diameter are required: one probe, called the 
ground probe, extends several centimeters below the liquid surface, while 
the second probe, called the surface probe, extends only 2 or 3 mm below 
the surface. The surface probe should be mounted no closer than 2 mm 
to the tank wall to prevent water drops from forming between the probe 
and the tank wall and giving erratic conductivity variations. 


BATTERY 



Figure 10. Liquid Height Gauge Circuit for Liquid-Slosh 
Quantitati ve Measurements 


The electrical conductance of the water between the probes is of the 
order of (10 kO)~\ and is largely independent of the distance between 
the probes. The liquid slosh measurements are based on the fact that this 
electrical conductance depends largely on the area of the surface probe in 
contact with the water. For a probe of uniform cross section, this area 
varies linearly with the depth of the surface probe ia the water. Then, for 
the circuit shown in Figure 10, the electric current flowing between the 
probes varies linearly with liquid height as does the resulting voltage 
across the 1000 signal resistor. 


219 



The voltage across the signal resistor thus reflects the changing height 
of the fluid at the surface probe due to liquid slosh. The surface motion 
should be restricted to ± 1 mm to ensure that the surface probe is always 
in the water, and to maintain a linear relation between voltage and surface 
height. Because of the small volume of each probe, several surface probes 
can be mounted at diflFerent places in the tank. These height gauges should 
be located so that different probes are sensitive to differeht modes of 
liquid slosh. 

R. A. Curtis of Comsat Laboratories has constructed and tested such 
a liquid height gauge. The apparatus gives satisfactory performance on 
the spinning test fixture described by Martin in Reference 8 (pp. 9-12). 
N. S. Lomas has used the fixture with a single height gauge and a tele- 
vision camera to look for resonances in a conispheric tank containing 
water at an 86-percent fill. A structural resonance limited the driving 
frequency ratio to values below 0.53, that is, to values far below the in- 
orbit resonance values. Liquid resonances associated with transverse and 
longitudinal oscillations were found for p = 0.349 and 0.407 respectively. 
These values are nearly one-half the resonance values of 0.72 and 0.80 
predicted for / = 0.86 by the straight lines fit to in-orbit data in Figure 9, 
The liquid was oscillating at the driving frequency. No firm conclusions 
about these newly observed resonances can be made because the height 
gauge showed many large-amplitude harmonics of the driving frequency 
whose presence is unexplained. 

Finally, experimental and theoretical consideration should be given to 
a possible swirl oscillation of the liquid, which is explained as follows. 
The body of the liquid, if treated as a rigid body sliding within a spherical 
tank, can oscillate with three degrees of freedom about three mutually 
orthogonal axes. Two of these axes are “parallel” to the equilibrium free 
liquid surface and correspond to the usual longitudinal and transverse 
oscillation of the liquid as a physical pendulum. The third axis is normal 
to the free liquid surface at the surface center and is the axis for swirl 
motion of the liquid. For a non-spinning tank on the ground, a displace- 
ment of the liquid about the swirl axis results in no restoring torque on 
the liquid; no swirl oscillation occurs in this case. On a spinning satellite, 
however, the equilibrium free liquid surface is a portion of a cylinder 
whose axis is the spin axis. A rigid-body displacement of the liquid about 
the swirl axis results in a restoring torque due to centrifi gal forces which 
attempt to bring the cylindrical surface back into alignment with the 
instantaneous s,p\n axis. The natural frequency of the resulting oscillation 
about the swirl axis may well correspond to an observed anomalous 
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resonance. This swirl oscillation would be driven by spacecraft nutation 
because the direction of the spin axis about the swirl axis oscillates with 
time. 
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The INTELSAT IV liquid slosh tests were a cooperative effort by many 
people at Comsat. The initial analyses that led to the tests were performed 
by L. B, Ricks and A. A. Satterlee 0/ Comsat and J. O. Salvatore of Hughes 
Aircraft Co. L. R Dest o/ CoMSAT helped reduce most of the nutation strip 
chart recordings and suggested many improvements to the text of this paper. 
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